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Abstract. Scattering of resonant radiation in a dense two-level medium is studied 
theoretically with account for local field effects and renormalization of the resonance 
frequency. Intrinsic optical bistability is viewed as switching between different spectral 
patterns of fluorescent light controlled by the incident field strength. Response spectra 
arc calculated analytically for the entire hysteresis loop of atomic excitation. The 
equations to describe the non-linear interaction of an atomic ensemble with light are 
derived from the Bogolubov-Born-Green-Kirkwood-Yvon hierarchy for reduced single 
particle density matrices of atoms and quantized field modes and their correlation 
operators. The spectral power of scattered light with separated coherent and incoherent 
constituents is obtained straightforwardly within the hierarchy. The formula obtained 
for emission spectra can be used to distinguish between possible mechanisms suggested 
to produce intrinsic bistability in experiments. 
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1. Introduction 

It has long been known that optical properties of dense atomic ensembles or complex 
materials may be very different from those exhibited by independent atoms. Some of 
these systems can produce two different output signals responding to the same input 
intensity from a driving laser. For the entire range of applied laser intensities the 
output function would form a hysteresis loop. When no external feedback is required 
for this phenomenon to occur it is recognized as intrinsic optical bistability (IOB). 
Experimentally IOB has been observed in the optical response from rare earth ions and 
ion pairs in glasses and crystallines (HEi[3j|IJ[5j[6j[7J[E]. However, its physical nature is 
still subject to debate. Most theoretical models suggest IOB originates from the fact that 
the transition frequency of a light emitter oja in a driven system is somehow renormalized 
and takes the effective value uja = — (,W . The latter is a linear function of inversion 
W where the renormalization constant Q follows from the coupling mechanism between 
the light emitter and its surrounding. Since W is uja dependent, it provides the system 
with an intrinsic feedback and nonlinearity and, therefore, produces bistability in the 
steady state solution. Most commonly such renormalization is introduced through near 
dipole-dipole interactions (NDD) between emitters [HI [TUJ, [TTJ. Alternatively, it can 
be achieved by considering model interactions for emitter pairs [121 EES EL coupling 
between a single emitter and the local vibration mode in its host medium [151 05] or in 
an interacting charge-transfer exciton system [17]. But, certainly, the IOB community 
has always been open to other ideas. 

Originally IOB was predicted by Bowden with co-workers for a dense ensemble 
of two-level atoms [I~8l [XQj [201 121] ■ Currently, in a simple semiclassical view to the 
problem, uja is a result of the famous Lorentz local-field correction (LFC). It states that 
in a dense homogeneous ensemble exposed to a laser light each atom is actually driven 
by an effective or local field strength 

47T 

E L = E M + — P, (1) 

where Em is the macroscopic field in the medium commonly approximated as the laser 
field and P is the macroscopic polarization. If one substitutes this condition to the 
optical Bloch equations this corrected field can be manipulated to give the effective 
ua = uja — (iW for which (l = (47r/3/i)iV|d| 2 is the Lorentz frequency or NDD 
interaction parameter. It depends on the number density of atoms N and the matrix 
element of the atomic transition dipole moment d. In the literature C,lW is often 
referred to as Bowden-Lorentz redshift [TO]. Thus it must be widely understood that the 
Bowden type excitation-dependent renormalization of the resonance frequency follows 
straightforwardly from the locally effective Rabi frequency Q = d • E^/Zi. Additionally, 
interpretations or enhancements of IOB involve cooperative processes [H [21 [22] and laser 
heating [6]. In particular, the work [22] pioneered the long series of papers studying local 
field effects in thin films. Various numerical models using nonlinear energy transfer 
processes in rare-earth-doped crystals have also been reported to show IOB behavior 
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[221 M 125] ■ 

Most IOB theories were built to describe either atomic excitation or radiation 
intensity hysteresis. However, not much attention has been paid to studying spectral 
properties of the IOB response. In the steady state limit the "classic" resonance 
fluorescence of a still two-level atom in vacuum is known to be a function of atomic 
excitation or inversion W. The line shape of the inelastic component of fluorescent light 
is determined by three frequencies: the spontaneous decay rate T, the detuning between 
the laser and the transition frequencies A, and the Rabi frequency Q (26J, [27J, EE]- It 
is apparent that as soon as the atom is placed in a host medium the frequencies take 
effective values T, A, and Q. Such metamorphose could be detected by an observer. 
We hypothesize that for an intrinsically bistable system with W hysteresis the spectral 
patterns would be evidence of which frequency is "responsible" . This is potentially a 
good technique to identify the adequate IOB model. 

Our work studies IOB as switching between fluorescence spectra and changing of 
spectral patterns with alternative mechanisms of bistability occurrence. As the first 
step we will study the original model of IOB based on the Lorentz local-field correction. 
We suggest a new method for calculations of emission spectra as a convenient and 
consistent addition to the conventional techniques [28]. The equations to describe the 
non-linear interaction of a two-level atom in the ensemble with the laser field and the 
properties of the scattered light are derived from the Bogolubov-Born-Green-Kirkwood- 
Yvon (BBGKY) hierarchy for reduced single particle density matrices of atoms and 
quantized field modes and their correlation operators. This method is a significant 
improvement to the atom-photon density operator formalism described in books [26l [29] 
and numerous papers. 

2. BBGKY and Master Equation 

We consider an ensemble of atomic particles and quantized field modes which we denote 
by subscript indices A and F respectively. The atoms are assumed to be "motionless" so 
that there is no interaction between them rather than via the field modes. The general 
structure of the Hamiltonian for such a system is given by 

H = Y,Ha + Y, H p + Y,Y, V ^ ( 2 ) 

A F A F 

where Ha and Hp are the energies of free particles and Vaf describes their binary 
interactions. In order to introduce our approach in a general manner we do not specify 
the explicit terms until we finalize the model. The evolution and properties of the 
system are found from the many-particle total density operator p(t) which obeys the 
von Neumann equation: 

iHj t p-[H,p] = 0. 

However, in our study we need to find the solutions for the reduced single-particle 
operators pA(t) and PF(t). The initial conditions for them can be defined using the 
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fact that if at t < no interaction is observed p(0) can be factorized to a product of 
all single-atom p A (0) and single-photon Pf(0) density operators. Further development 
of the single particle operators as constituents of the "unfavorable" system and the 
steady state mode are to be revealed. This goal can be approached if we substitute the 
von Neumann equation with the BBGKY hierarchy for reduced density and correlation 
operators [301 Ell E2] ■ Some issues of building such an hierarchy for atom-photon systems 
and proper decoupling of equations were discussed in [33j EI]- In this paper we will skip 
the cumbersome derivations and only point out the main steps to the analyzable system 
of equations. First of all, let us note that in order to treat the problem properly one 
needs to consider the actual number of particles and avoid using the thermodynamic 
limit where the problem is treated in terms of constant densities of material particles. 
We formally assume all atoms to be individuals and turn to their bulk concentration 
only when performing spatial integration. Correspondingly, the reduced matrices are 
obtained from p by taking a partial trace, i.e. 

P{S} = Tr {S '^s}P, Tr {s} p{s} = 1. (3) 

Here, {S} and {S'} denote non-overlapping collections of species such that {S' + S} is 
the system described by p. The reduced multi-particle matrices are substituted with 
their cluster expansions [301 El] : 

Paf = PaPf + 9af, (4) 

Paff' = PaPfPf' + 9afPf' + 9af'Pf + Qff'Pa + 9aff>, 

etc. 

From this definition it follows that in order to meet fl3]) the trace or even a partial trace 
of the correlation operator g^s} should be strictly zero: 

Tr {s}9{s} = 0. (5) 

With all these starting points the von Neumann equation could be equivalently replaced 
with an infinite series of coupled kinetic equations. The first two types of equations are 
precisely 

d_ 

"Jt 1 

d 
'~dt' 



ih di PA ~ ^ Ha,Pa ^ = Trp [Vaf > 9af] > (6) 

F 

ih Ti pF ~ pF ^ = Tr A v AF, 9af]- (7) 



A 

Here, the l.h.s. terms constitute the von Neumann equations for single-particle operators 
coupled on the r.h.s. to each other and the remaining system through the atom-field 
correlations found from 

ifi—g AF - [H A + H F , g AF ] = [V A f, PaPf] (8) 

+ ^2 Tr F"l V AF',9AFF']- 
F'^F 
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This equation is written in the limit of the generalized second Born approximation 
[3Tj 134"] for which the contributions from other two-particle correlations are neglected 
and only partial account for the three-particle correlations is taken: 

ih^g AFF , - [H A + H F + Hp,, g A FF>] (9) 
= [Vaf 1 , Pf'9af]- 

The primes at the particle-reference indices are used to distinguish different particles of 
the same sort. The higher order correlations are beyond the scope of the current analysis. 
Equations (JH])-([9]) are closed and energy conservative. The single particles here are also 
coupled to the remaining system via the semi-averaged interaction operators or Hartree 
(mean field) terms in the effective Hamiltonians: 

Ha = Ha + / ](Vaf)f, 

F 

H f = H f + Y,(Vaf)a, (10) 

A 

Vaf = Vaf — (Vaf)f — (Vaf) a- 

The averages must be understood as (X)s = Trs(Xp). Note that the interaction 
operator ( fit)]) corrected with the Hartree terms is the one that enables the correlation 
operators be strictly zero when traced over any correlating particle as stated in (JS]). 
Now we can specify the light emitters as two-level atoms and for free particles write 

H A = huj A or z A , Hp = hu k al s a Ks , (11) 
The atom-photon interactions will be treated in the electric-dipole approximation 

V A F = -dA-E K s(r A ), (12) 
for which the dipole moment operator of an atom localized at some point is 

d A = d A a A + d A a\ (13) 
and a single mode of the electric field at r is given by 



E k .„(r) = ^/^^e Ks d Ks e^ + H.c. (14) 



Here, a A , a A are the population inversion, raising and lowering operators for atoms with 
transition frequency to a and is its dipole transition matrix element. Then a\ iS and 
ak, s are the creation and annihilation operators that refer to a field mode with frequency 
Uk, wave vector k, k = Wfe/c, and polarization e^ s . V denotes the quantization volume. 
The action of the classic external field is reflected by setting the initial condition for a 
number of selected modes: 

p F (0) = \a k:S )(a Ks \, (15) 

where |ak, s ) indicates a coherent state of the field. We will further use the space-time 
commutation relation 

[E(r, t), E(r', t')} = -iKG{r -r',t- t') (16) 
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for the electric field 

V f 

E(r,t) = 5> k , s (r,t), E^7^M k E' 

where G(r — r', t — t') is the Green's tensor following from the dyadic product of the 
field operators [35j 121] • At the same time the dyads can be written as 

E(r, t)E(r', if) - I(r, r'; t, t') = -ihG + {r -r',t- t'), 

E(r',t')E(r,t) - i(r,r';t,t') = zftG _ (r - r', * - if), (17) 

for which we introduced the field correlation tensor 

i(r,r';M') =: E(r, t)F>(r',t'): 

where :: denote normal ordering and are the retarded and advanced Green's tensors 
respectively. The important averages we use below are 

I(r,r'-t,t') = (I{r,r'- t,t')), 
I inc {r,r';t,t') = ((I(r, r'; t, t'))) 

= (:E(r,t)E(r' ! 0:)-(E(r,t))(E(r',t')). (18) 

Equation ([6]) can be shown to give the master equation if one eliminates its field operators 
using the standard techniques of quantum optics. First, it is beneficial to transform to 
the interaction picture for which 



p = exp 



cxp 



In this picture the remaining Hamiltonian operators (TLB"]) and (|14p describe interactions 
and acquire the explicit exponential time dependencies ±iu A t and iiu^t. Now we can 
substitute the formal integrals of ((Tj) and ([9]) over time t' < t into ([6]). Then we assume 
p(t') ~ p(t). Retaining the terms which contain |d^| and |d^| 2 only and making the use 
of (}TE]) - (JT8]) with r' — > r A we come to 

^Pa = | d A (t) ■ E L (r A ,t),p A 

+ i jf <Hf[d A (t), ■ G + (0, t-t')- d A (t')p A ] 

+ % - y* dt'[d A (t),-G-(o,t-t')-p A d A (t')} 
i r* 

+ -J dt'[d A {t),[-I inc {r A ,r A ;t,t')-d A {t'),p A }}, (20) 



where 



El — E 



+ N f dv A , [ dt'G(v A -r A ,,t-t')-(d AI (t')), (21) 
Jv Jo 
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The first commutator on the r.h.s. of (I20p describes interaction of atom A with the 
local field (l2Tj) in which E is the external field. It originated properly in accordance 
with the condition defined in ( IT5l) . The integrals in ( [2T]) describe contributions from 
the other atoms in the ensemble to the resulting effective field acting at position r^. 
Here we used the continuous medium approximation and replaced the summation over 
atoms with a proper integration over the sample volume V. The prime reflects the 
property of BBGKY which excludes atom A from the mean field ffTUl) and eliminates 
its contribution to ( j2~Tj) . We will apply the conventional approach to such integration 
where one should mind an "empty" sphere around atom A. 

The next two integrals in ( 1201) describe atomic relaxation due to spontaneous 
emission. The last term makes account for atomic transitions induced by incoherent 
photons. However the latter processes are to be neglected in this paper. In order to 
make a comparison with the "classic" works on IOB we restrict our model to the case 
when there is no reabsorption of scattered light. 

The equation for the local field does not make an account for the complete Green's 
tensor which has the delta-function peculiarity at zero point [33J ESI [37] because this 
region has been removed from the integration volume. However, if we rewrite ( f2~Ti) with 
the full tensor and separate integration over volumes V and 5V = V — V the local field 
becomes easier to handle with for homogenous and isotropic media: 

El = E 

/"* f An 
+ dt' dr A ,G{r A -r A ,,t-t')--p{t') + —P, (22) 

Jo Jv 6 
where we introduced macroscopic polarization of the sample P = N(d A >). Equation (T2"2"j) 
is exactly the well-known Lorentz local field condition ([1}) which describes the difference 
between the macroscopic field in a medium and the field acting on each atom therein. 
This is the first result of this paper which demonstrates that the master equation making 
an account for the local-field correction can be self-consistently obtained from BBGKY. 

3. Bistability and emission spectra 

The master equation obtained in section [2] can now be rewritten in a more transparent 
form. Without its last term (120 p becomes independent and may be easily transformed 
to the IOB equation using the Born-Markov and the rotating wave approximations. 
Changing back to matrix p in ( 1T9|) and operating in the frame rotating with the laser 
frequency ul we get 

^p A = - ?[Acr z , p A ] + i[a + Cl* + p A ) 

- ±T {[&+,&- p A ]- [*-, PAO-+]}, (23) 

where A is the detuning, Cl is the effective Rabi frequency and T is the spontaneous 
decay rate. Since all operators and frequencies are associated with a single arbitrary 
atom we may now drop the unnecessary subscripts to simplify the notation. Given 
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the model described above these parameters are exactly the following: A = uja — ^l, 
fl = d • E L /h, and r = 2/h d • hnG(0,u A ) ■ d = (4/3ft)|d| 2 fc 3 . The contribution of the 
Lamb shift is neglected. 

This form of the master equation has been very well studied in the steady state 
limit and is known to describe IOB (see section [TJ. When ( 1231) is written as the set 
of equations for matrix elements and resolved for the population difference function 
W = —2{a z ) it produces the third-order algebraic equation: 



Coefficient ( in our case is precisely the NDD parameter that follows from ( 122]) . 
Consequently, we are to write Q = Q + C(^" + )) where Q = d ■ Eo/h and C = (l = 
(47r/3/i)iV|d| 2 = nNT/k 3 . Three real roots of ( I24p found for various laser field strengths 
(~ Q) form the hysteresis loop of atomic excitation p 2 2 = (1 — W)/2 as shown in 
figure [TJ Similarly, it is also the function which describes IOB in terms of the total or 
frequency- integrated power of fluorescent light as A mt ~ Tp 2 2- The coexisting steady 
state solutions are notated as p l $° and p 2 f an d refer to the lower and the upper curves 
respectively 

Equation (121"]) can as well be derived from ( 123]) in which A — > A = A — (W and 
Cl = Q. As discussed earlier it could be, for example, either the case suggested in [T5l[T6] 
or in [IT]. In this paper we will not reproduce those models using BBGKY formalism 
as it would require chains of cumbersome transformations similar to those outlined in 
section [2j With no new features revealed it could be reasonable to introduce the effective 
detuning phenomenologically. In this approach we assume C = Cm = mY, where m is 
an arbitrary factor. In any case with either effective detuning or Rabi frequency we use 
( l24l) with ( = (l = Cm to find the source of fluorescent light. 

Calculation of fluorescent spectra requires the use of (J7J) and (jHJ). In [331 El] the 
spectrum was defined as irreversible loss of radiation in a nonconservative system. It was 
shown that such treatment was necessary to calculate the time-dependent or transient 
mode spectra correctly. The spectra relevant to the steady-state hysteresis do not 
demand this kind of special care and can be evaluated within the set of equations (jSJ)-® 
that conserve energy. However in this paper we will follow the approach of [33| 131] to 
demonstrate succession of the research techniques applied for related problems. 

Let us note first that elimination of three-particle correlation qaff 1 in (02) produces 
operators similar to those found in r.h.s. of (|20|) . Dropping the terms dependent on the 
number of photons we get the atomic contribution to the correlation function associated 



c 3 W 3 + c 2 W 2 + c x W + c 
where the polynomial coefficients are: 



0. 



(24) 



c 3 = C 2 , 

c 2 = -C(C + 2A), 
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Figure 1. Excitation hysteresis as function of the laser field strength for A = 3r and 
C = 50I\ Solid lines indicate the steady state solutions, dotted line is the unstable 
branch. The dash-dot line describes the saturation curve of a single independent two- 
level atom. The laser field strength is expressed in terms of the Rabi frequency f2 in T 
units, P22 - population density of exited atomic level. 

with the Lorentz local field and spontaneous emission. Besides, we will supply © and 
()9]) with the field damping operator: 

£ kjS X = -%/2{[aj c s , a kjS X] - [a Ks , Xa^J}, (25) 

where 7] k is the loss rate at zero temperature. This damping could be introduced 
accurately within BBGKY if the system is completed with the field reservoir. However, 
since it is a well known result we simply write it in. Later on in our calculations the 
damping rate is considered to be rjk —> 0. 

Meanwhile we have a complete set of equations to study spectral properties of 
resonance fluorescence. The power density of a field mode is a function of 

A' k = huj k ^-(al s a Ks ) = hu k Tr F ^al s a k , s ^p^|- (26) 

The time derivative of pp is defined by (J7J) with attachment of (125]) . Since we are to 
count photons leaving the system we need to consider the nonconservative term only. 
The latter is nothing but ( |25i) . It is also necessary to write the opposite sign to become 
an "observer" . In such a way we come to the function that describes the power of the 
emission signal absorbed by an idealized detector: 

A k = -fw k Tr F {-ial s a k!S C kiS p F } = hu k rj k (a\, s a Ks ) . (27) 
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It must be understood that the photon counter is characterized by a continuity of rj k . 
Now our task is to find rj k (a/ k s ak,s) m the steady state. The proper equation for the 
number of photons follows straightforwardly from (J7J): 

A 

+ J2 Tr ^ + 9k) + H.c. (28) 

where g k = ^(^^fA^af)- The atom-field coupling is represented here by 
the coupling constant £fc( r ) = \f2-Ku k /fiN{d ■ e^ s )e tkr . In the steady state with 
d(d k s ak, s ) / 'dt = evaluation of A k becomes obvious. Functions {d^ s ) and {d ks ) are 
eliminated naturally using ([7J and ( 1251) transformed to the form 

j t « s ) = - *f )<a k , s > + J]CM(r). (29) 

Function and its conjugate are found from the equations obtained from and fl25|) : 
+ i^fc + *[A«t«, </fc ] - i[cr + tt* + gl\ 

- \{[^^-9t\ ~ lt-,9t* + \} = M\o-p A - (a-)p A ), (30) 

where we introduced the scale frequency v k = w^-wl and applied rj k — > 0. In derivation 
of these equations we always assumed that all the approximations used to obtain the 
master equation (1231) were also valid. In the steady state for N light emitters with no 
account for the spatial part we get 

Vk(al s a kiS ) = 2iV 2 ^|p 12 | 2 ^ + N{(gt) 12 + H.c). (31) 

Finally using 

lim - ^ = 6(x) 

to substitute (13T1) into ( [271) the expression for the fluorescence spectrum becomes 

A k = N 2 F k \p 12 \ 2 S(u k ) + 2iV^ fe Re{(<?+) 12 } (32) 

where the off-diagonal matrix elements pi 2 = (<x + ) and (g k )i2 = TrA(& + g~l) are yet 
to be found from ( 1231) and ( 1301) respectively with the time derivatives put to zero. F k 
is a power density dimension factor. It is important to note that ( 130]) is linear and 
can be solved explicitly as shown further. When Q = Q the final expression becomes 
the famous formula for the fluorescence spectrum with coherent and incoherent parts. 
However, unlike the method described in [261 129] we performed no artificial steps to 
separate the two parts. Therefore, the technique we present is convenient for studying 
spectral properties of systems exhibiting nonlinear response to external light for which 
analytical solution of ( 1231) is hard or impossible to find. 
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Figure 2. Emission spectrum at point 1 (lower branch) for A = 31", 51 = ffi ks 15. 6I\ 
The solid line represents emission spectrum for £l = 50r, £ m = 0, Dotted line - for 
Cm = 50r, Cl = 0. The gray shape describes the emission spectrum for an independent 
two-level atom. A® is maximum value of the emission power function for a single 
independent atom in the saturation limit. 



4. Comparison of spectral patterns 

As we have already pointed out the equations obtained in the previous section are 
suitable to describe two mechanisms of bistability, i.e., IOB produced by either Q or 
A. The latter model will give (130]) effective detuning A = A — ( m W and Q = Q. 
In this section we will show the differences in spectral patterns of inelastic fluorescent 
emission expected for Q and A based IOB models given (l = Cm and excitation curves 
are identical as plotted in figure [TJ We can perform a formal generalization of (130]) by 
writing Q and A and deciding which is to represent an effective value later in the final 
expression. From (132]) it follows that we need the function of Re j (5^)12 j which can be 
found from (130]) . It can be expanded to the system of four equations for four matrix 
elements (9^)11, (5^)12, {9^)21, (6^)22- However, we can reduce the number of equations 
to three using the general property of a correlation operator (J5j) which in this case gives 
(9k)u + (5^)22 = 0. Consequently, the set of equations can be represented in the matrix 
form as follows: 



M x g = q, 



(33) 
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where M is the matrix: 

/ ivk + r 



M 



in* 



2iQ i{y k - A) + \Y 



-ifl 







(34) 



V -^n* o i(u k + a) + |r / 

and g and q are vectors. The first represents the correlation function in terms of its 
matrix elements 

T 



9 



(0fc)n. (9t)i2, (gthi 



(35) 



Vector q = \(,k\ 2 s — |£fc| 2 /?2iPA has two constituents where s represents the source of 
spontaneously emitted photons 

s = ( P21, P22, ) (36) 
and vector eliminates the Rayleigh singularity 

Pa — ( Pu, pn, P21 ) ■ (37) 



Solution of fl3"3"|) gives 

(9k) 12 = 



~2p2 2 l6c| 2 



v\ - T 2 - 2|fi| 2 - HYv k 
detM 



(38) 



In order to obtain the desired expression for Rej^^)^ j let us multiply the numerator 
and denominator by the complex conjugate term det M* and take the real part of the 
numerator: 



Re{(^ + ) 12 }=2p 2 2 |e fc | 2 r c 



a 



^ + 64^ + 63^ + 60 
The corresponding polynomial coefficients are: 

1 



(39) 



a 
b 4 . 



2\Q\ 2 + A 2 + -r 2 
4 

= 2|fi| 2 + r 2 



-8101 



2A 2 + -r 2 

2 



16|fi| 2 + 2|fi| 2 (4A 2 + T 2 ) + A 4 

1M21QI- + A 2 + -r 2 " 

4 



3 r 2 A 2 + _9_p4 

2 16 



Equation (I3"9l with proper substitutes for either Q or A is the main results of our paper. 
It determines the shapes of the inelastic emission spectra in a wide range of parameters. 
Below we present several plots of the spectrum profiles corresponding to four important 
points on the hysteresis curve in figure [TJ Points 1 and 2 give the atomic excitations 
around the threshold value of Q = that refers to the switch from the lower to the 
upper branch of the hysteresis loop. For a single two-level atom in the steady state 
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one cannot expect the excitation to be as low as indicated by point 1 in the limit of a 
strong field. The threshold point for the lower curve is sufficiently far in the saturation 
region. Consequently, the spectral function ( 139]) for p 22 = p^ffi) is notably different 
from conventional triplets [27J|28]. Figured demonstrates the plots for Q and A variants 
of IOB. One can see that the line splittings and peak values are sharply "outlined" by 
the respective models. The distance between the satellites and the central peak is over 
10 times greater in terms of T units for IOB response conditioned by effective A as 
compared to the Lorentz model. Moreover at this point the linear response spectrum is 
sharply distinctive as it is nothing but the "classic" Mollow's profile. 

1.0 I 1 1 

o ^ 

0.8 - 

< 



O) 0.6 - 

'(f) 




Frequency v k IT 



Figure 3. Emission spectrum at point 2 (upper branch) for A = 3I\ Q, = m 15. 6I\ 
Solid line represents emission spectrum for £l = 50r, ( m = 0, dotted line - for 
Cm = 50r, (l = 0. 

For the upper branch at point 2 we have p 2 2 = P 22 (^) an d the emission spectra 
plots are shown in figure [3j Since here W ~ 0, the contribution from the effective 
detuning is weak, i.e. A « A, and Q » A,T. Obviously, the respective triplet 
coincides with the classic shape. In its turn the Lorentz IOB spectrum still reflects some 
difference between Q and Cl. However, it will gradually vanish with further increase of 
Q. This can be very well seen from the strict relation 

\tt\ 2 = \tt\ 2 A2 + r2//4 (40) 

that follows from (T23|) . As one can see in (|39|) the line shape depends on \Q\ 2 and the 
relation to the spectral properties of a single atom is determined by (HO]) . 

As we reduce the strength of excitation we move to point 3 on the upper steady 
state curve. It corresponds to the weak threshold field Q^, high excitation and small W. 
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For these parameters the emission spectra are shown in figure HI Here we find smaller 
splitting between spectral lines as compared to the single atom picture. Opposite to 
point 1 the A bistability triplet has a narrower range in comparison with Q IOB. 
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Figure 4. Emission spectrum at point 3 (upper branch) for A = 3r, f2 = fir « 1.6r. 
Solid line represents emission spectrum for £l = 50r, £ m = 0, dotted line - for 

C m = 5or, a = o. 

Further reduction of the pump strength will make the system transit to point 4 at 
which W ~ 1. Since Ak is proportional to {p l ^"{^)} 2 the signal intensity should be 
very small and disadvantageous for registering the line shapes. 

Let us consider the denominator in equation ( )39|) for the emission spectrum: 

4 + h»t + b 2 vl + b (41) 

It can be written as a product of the peak roots and two positive terms in the following 
form: 

' (42) 



where 



and 



u„ - \/4|fi| 2 + A 2 -^r 2 (43) 

7 6 = T 2 ^2|fi| 2 + A 2 + (44) 

From (|42j) the triplet structure of fluorescence becomes clear. The central peak at v^. = 
is symmetrically surrounded by the satellites with maximums at v k = ±z/ p . It follows 
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Figure 5. Positions of emission spectra peaks for the lower and upper branches of 
the hysteresis loop. Solid line represents peaks position due to local field correction 
Cl = 50r, Cm = 0, dotted line - due to effective detuning ( m = 50T, (l — 0. Gray solid 
lines describe positions of spectrum peaks of an single independent two-level atom with 
a laser field applied in the vacuum. 



straightforwardly from the fact that if (Vfc±z/ p ) 2 = function ( 1391) reaches its maximum 
value. In addition to the conventional analysis of fluorescence these expressions contain 
the effective or renormalized frequencies. Equation ( 1431) gives the exact position of 
the side peaks for both models of IOB. One may write the proper substitute for \Q\ 2 
or A 2 to introduce the respective mechanism and obtain a picture similar to the one 
shown in figure [5j It is demonstrated here that the frequency shifts of the satellites 
in a strong fluorescence signal may vividly characterize which mechanism is likely to 
produce a saturation curve with an IOB hysteresis loop. Moreover, ( H3|) is also suitable 
for studying joint contribution of effective parameters. However, this issue as well as 
calculation of p 2 2 to find a complex IOB picture is out of discussion in this paper. 
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Finally, we may note that the Bowden-Lorentz type of IOB is best detected at point 3 
or, alternatively, on the lower branch as narrowing of the spectral range. At the same 
time the real renormalization of the atomic transition frequency could be very well seen 
when uprising to the threshold point 1. 

5. Summary 

We have performed a purely analytical study of IOB as switching between the patterns of 
fluorescence spectra. The local field induced hysteresis was shown to originate naturally 
for the model of a dense collection of like emitters when described thoroughly within 
the BBGKY hierarchy. The steady state spectra were found self-consistently within the 
developed formalism using the atom-field correlation matrix. We demonstrated that 
in terms of spectral analysis the contribution from the Lorentz local field cannot be 
described via introduction of an effective excitation depended detuning. In other words 
no evidence of actual renormalization of the atomic transition frequency can follow 
from observing light scattering in a monocomponent system. Within this concept the 
incoherent Mollow triplets would reflect abnormally small effective pumping of light 
emitters along the lower branch of the hysteresis loop. At the same time, a significant 
increase of the field induced splitting must be expected for the high excitation state. 
The latter though should gradually disappear in the saturation limit. 

In order to make a comparison with the models suggesting renormalization of the 
transition frequency due to interactions of the atom with its complex environment 
we introduced this property phenomenologically. Since the IOB equation for atomic 
excitation is valid for both concepts the system parameters were chosen to give identical 
hysteresis loops. Analysis of the spectral profiles for this case revealed substantial 
differences in the character of line splitting at weak and moderate excitation. The 
results of our comparison are supported by analytical equations fl39l) - PU|) and f )43|) . 

This analysis was carried out to emphasize the importance of spectroscopic data 
in studying optical response of complex materials. The problem could be described 
theoretically with the BBGKY formalism adjusted for atom-photon interactions. It is 
very advantageous to use this approach because the theory could be created from the 
first principles avoiding common phenomenology. The BBGKY method applied in this 
paper gave us a self-consistent set of equations with remarkable features which revealed 
some optical properties of the atomic system and the scattered field. It was due to the 
BBGKY concept and the general property of correlation function that we succeeded in 
reducing the number of equations to obtain the fluorescence spectrum function. The 
Rayleigh elastic constituent of the scattered light was automatically separated from 
the broad inelastic components which simplified our calculations greatly as compared 
to other spectroscopic techniques using the density matrix. The formula obtained for 
emission spectra can be used to distinguish among possible contributions from effective 
atomic frequencies in a medium showing IOB in experiments. 
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